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ABSTRACT 

If  f^^(a)  (a  =  a,b,  i  =  0,1,..., k  -  1)  are  given,  then  we  get  a  class 
of  the  Hermite  approximation  operator  Qf  =  F  satisfying  F^(a)  =  f^(a), 
where  F  is  the  many-knot  spline  function  whose  knots  are  at  points 
yi  :  3  =  y0  <  Y1  <  ***  <  Yk-1  =  b/  3nd  F  e  Pk  °n  tyi_i'yil*  The  operator 

k- 1  /  •  \  j  j \ 

is  of  the  form  Qf  :=  £  [f  (a)$  +  f'  ;(b)^.l.  We  give  an  explicit 

i=0  1 

representation  of  <J> ^  and  in  terms  of  B-splines  N^  We  show  that 

Q  reproduces  appropriate  classes  of  polynomials. 
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Key  Words:  Hermite  interpolation,  splines,  local  explicit,  many-knot 
spline  function. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  deals  with  Hermite  interpolation  on  the  interval  [a,b]  using 
many-knot  splines.  The  contribution  of  this  paper  is  to  find  a  many-knot 
spline  function  F  of  degree  k  -  1  whose  knots  are  at  points  y^: 

a  =  yQ  <  y,  <  ...  <  y^,  =  b  . 


When  conditions  are  given  on  the  ends  of  [a,b]  :  f^^(a),  <*  =  a,b, 

k-1 

i  =  0,...,k  -  1,  then  F  =  Qf  =  £  tf  1  (a)f  +  f  1  ( b )  .  ]  and  F  e  P  on 

i=0 

[y  ^  ,y .  ]  i  =  1 ,2 , . . . ,k  -  1 .  The  explicit  representations  of  basic 
functions  <J>^,  which  have  properties 

<f>(l)(a)  =  6..,  <j,(l)(b)  =  0,  i|»(l)(a)  =  0,  i|>U)(b)  =  6.. 
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for  all  i,j  =  0,1, ...,k  -  1 


are  given  in  terms  of  B-splines  We  also  prove  that  this  approximation 

operator  Q  reproduces  appropriate  classes  of  polynomials. 

Since  the  degree  of  the  many-knot  splines  used  here  is  lower  than  that  of 
ordinary  Hermite  interpolation  and  the  knots  of  the  splines  can  be  chosen,  it 
would  be  useful  for  some  problems,  for  example,  in  Computer  Aided  Geometric 
Design  (CAGD).  r  - 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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\  ■  ‘ 

1.  \INTRODUCTIO * 

> — i  — . . 

r -r  ,  , 

Soae  authors  considered  operators  of  the  form  Qf  ■  l  A  fN  where  tN  .  t  is  a 

Z,,  1  1  f  K  1  $  K 

sequence  of  B-splines  and  {A^}  is  a  sequence  of  linear  functionals.  The  variation 
diminishing  method  of  Schoenberg  ([9),  (5),  161),  the  quasi-interpolant  of  de  Boor  and  Fix 
are  well-known.  Such  approximation  schemes  have  several  important  advantages  over  spline 
interpolation.  They  can  be  constructed  directly  without  matrix  inversion,  local  error 
bounds  can  be  obtained  naturally.  Qi  considered  so-called  many-knot  splines  which  have 


many  more  knots  than  degrees  of  freedom  and  constructed  the  cardinal  spline 

»  -  1  ,  , 

Qf  »  >  f(x,)q,  ,  where  q,  v  is  made  up  of  B-splines  on  a  uniform  partition,  has  sna 

support  and  satisfies  q,  ,<x , )  «  d  •  ‘  Such  an  approximation  operator  reproduces 
l,k  j  ij 

appropriate  classes  of  polynomials^ 

The  purpose  of  this  paper  is  to  construct  a  class  of  many-knot  explicit  local 
polynomial  spline  approximation  operators  for  Hermite  interpolation  of  real-valued 
functions  defined  on  some  interval  [a,b]. 
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Let  be  the  set  of  polynomials  of  degree  less  than  k,  and  let 


*0  <  *1  < 


<  *k- 


k-1 


b  . 


(1.0) 


We  define 


'"’'-.v"*'1"'' . “2>' 


is  the  familiar  class  of  polynomial  splines  of  order  k  with  knots  at  the  points 
(i  *  0,1, ... ,k  —  2). 
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Let  F  be  a  linear  space  of  real  valued  functions  on  [a,b),  and  suppose  F 
contains  the  class  of  polynonials  P  .  Given  f  e  F,  we  construct  an  approximation 
F ( ♦ )  ■  flf ( • )  such  that 

f(t)(a)  -  f(t,(a).  F(1)(b)  -  f<4>{b)  t  -  0, 1 , . . . ,k  -  1.  (1.1) 

In  other  words,  set 

k-1  ...  k-1  ... 

Qf  s-  l  f‘3,(a)4  (x)  +  l  f'3'(b)4  (x)  ,  (1.2) 

j-0  3  j-0  3 

suppose  4^  satisfying 

*ji><a)  "  Sij  '  “  0  (1.3) 

4*i,(a>  -  0  ,  (i^tb)  -  6  (1.4) 

i, j  ■  0,1, ...,k  •  2  • 

If  4^  and  4^  ere  chosen  in  *2k-2'  then  thi*  problem  above  has  been  considered 
(see,  for  instance,  [1],  [31,  (4]),  and  in  this  case  F  e  pJk_2  on  [“»*>)• 

He  will  find  a  many-knot  spline  P  e  satisfying  (1.1).  Such  many-knot  cardinal 
splines  (4^)  and  (4^)  are  of  degree  less  than  k,  therefore  r  is  also  of  degree  less 
than  k.  He  present  4^  and  4^  as  explicit  representations. 

This  paper  proves  that  the  many-knot  spline  Heraite  approximation  operator  Q 
reproduces  appropriate  classes  of  polynomials  on  [a,b]. 
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In  order  to  consider  the  general  case,  denote 

I 


{0,1,...,n} 

VX>  *  °j.MxU  '  *  *  tXi'Xi+l’ 

j  »ei  ■“ 


k- 1 


1  e  *k-2  '  j  e  Xk-2 


(the  partition  is  0  »  xQ  <  x1  <  <  . . .  <  xk_^  “  1), 


and 


♦‘1,(xi  -  0)  »  +0),  i!  Ik_2\{0).  t  e  Ik_2  , 

♦5A><0)  -  «iJ  ,  ♦^i)(i)  -  0  ,  i,j  e  Ik_2  . 

Since  we  have  k(k  -  1)  unknown  coefficients  a  with  k(k  -  1)  conditions,  so  it 

j » 1* 

seems  possible  to  find  a,  •  But,  it  is  difficult  to  get  the  explicit  representations  for 

j/U 

Below  we  will  directly  present  the  explicit  formulas  for  and  i|>^. 

Here  are  the  notations  used  in  our  discussion. 

Let  X  i»  (x^)  be  a  nondecreasing  sequence.  The  i-th  B-spline  of  order  k  for  the 
knot  sequence  (x^)  is  denoted  by 

■i.k00 !“  <xitk  •  xi,txi . w‘-»r 

for  all  x  e  R,  where  the  symbol  (x^, . . . ,x  )  denotes  the  k-th  oder  divided-difference 
functional 


*VVaa . Vi’  !“  i  “v  %  •••  av 

(V1, — ,vy)  V1  v2  U 


v,  e  i 


n-1  '  Wi  *  V1  *  J» 


r(H)  _  ,  .  ,(K  -  1 ) 

«i  ■y"M.1<xi+1»xi+1 . xi+k.,>/lu  - 


j-  symQ(...)  :»  1  • 


Proa  (1.0),  we  define 


*i-  1 


‘i  ' 


fi  “s  Xk-1+i  '  for  1  e  Xk-1 


(2.2) 
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Thus  wo  get  a  partition  on  [-1,1]  from  [0,1] s 

-1  -  x0  <  x,  <  ...  <  xk_2  <  xk_t  -  0  <  xk  <  xk+t  <  ...  <  x2(k_1(  -  1  . 

We  construct  the  following  functions  on  [0,1]  as  a  special  kind  of  combination 
B-splines 

v«  -  37 . • 

k-2 


for  x  e  [o,i]  ,  j  e  ik_2  . 

Theorem  1.  The  functions  ♦^(x)  defined  in  (2.4)  satisfy 

*  4tj  . 

*^’(x)  «  0  for  |x|  >  1  ,  t,J  e  lk_2  . 

Proof.  If  i  e  I  and  |x|  >  1,  then  n!**(x)  »  0,  therefore  ♦***(x)  «  0 

1|K  j 

t.j  e  lk_2  and  |x|  >1.  If  i  5  I^,  then  N^fO)  -  0  since 

ik2  -  {ill  e  {...,-2.-1,0,l,2,...),Nifk(0)  *  0}  . 

By  Marsden's  Identity*6!,  for  x  e  [0,1] 


_U-1 


I 

i«I 


tiM,*,i,k(X)  '  M  ”  1'2' 


,k  . 


2k-3 


Thus 


x-0 


3  'x-0  5  ielX-2  ^  i,X 

hm  +2TKl+,\k(x,](“| 

3  iei,.  _  i-k-1  1  1 


k-2 


Let 


-  jj-  (x3)(t,|  -  ,  for  t,j  e  I 


♦j (x)  j-  ^  (x  -  1) 


k-2  * 


(2.3) 

of 

(2.4) 

(2.5) 

(2.6) 

for  all 

(2.7) 
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In  uniform  case 


♦o  =  N0,4  +  N1,4  +  N2,4  ' 


>1  ’  -TN0,4  +TN2,4  ' 


x  e  [o,i] 


*2  =  54  N0,4  ”  54  N1,4  +  54  N2,4 


3,  THE  OPERATOR  Q  REPRODUCES  APPROPRIATE  CLASSES  OF  POLYNOMIALS 

Using  the  functions  $ .  and  #  we  have  the  following  approximation  operator 

QfC)  :=  l  [f‘ j)(0)*.  +  f<j)(1)*.]  (•)  , 

iei  3  3 

J  k-2 

Q  defines  a  linear  operator  mapping  F  into  S^. 

Theorem  2.  Qg  -  g  for  all  g  e  P^. 

Proof.  Let 

span  (N)  span(Ni  k;  i  €  I2k_3>  > 
span (♦,’]')  :=  span  ( ♦ ;  j  e  Ik_2>  < 

S  :=  {g  :  Qg  =  g}  . 

Then  both  span(N)  and  span($,i(i)  are  linear  subspaces  of  F  on  [0,1]  of  dimension 
2k  -  2. 


By  (2.4)  and  (2.8)  we  have 


Since 


span($,i|i)  C  span  (N)  . 


Obviously 


i.e. 


dim( span (4 , V) )  =  dim( span(N ) )  =  2k  -  2  , 
span(^,\)i)  =  span(N) 

£  span(N) 

C  span(^,+) 
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S  »  span  (♦,!)>) 


(3.1) 


Now  it  ia  sufficient  to  prove  that 


It  follows  from  the  definition  of  the  set  S  and  the  operator  Q  that 

S  C  span( .  (3.2) 

On  the  other  hand,  Theorem  1  implies  that  we  have  Qf  =  f  for  any  f  e  spanl#,^). 

Hence 

span!*,*)  C  S  .  (3.3) 

(3.2)  and  (3.3)  mean  that  (3.1)  is  valid. 
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